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Abstract 

We have investigated two flavor neutrino oscillations in medium 
with periodic step electron number density profile. An approximate 
analytical solution have been found when the length of the density 
fluctuation is smaller then the neutrino oscillation length. 

Introduction 

Neutrino oscillation mechanism gives an explanation for decreasing flux 
intensity of electronic neutrinos coming from the Sun. The problem was 
predicted in p] and was solved via averaged vacuum oscillations with large 
mixing angles in [2]. In [3], [I], [5] solution for neutrino oscillations in matter 
with constant electron number density were obtained. An adiabatic neutrino 
oscillations were investigated in [6], [7]. An exact analytic expression for the 
probability amplitudes describing the oscillations of two (flavor) neutrinos in 
matter with exponentially varying density was obtained in [8] . This theoret- 
ical picture of such potential is in good agreement with the standard solar 
model according to which the matter density in the sun decreases exponen- 
tially with the distance from the center of the sun, except for the regions 
located close to the center and to the surface. But such a picture doesn't 
take into account the unpredictable variations of electronic gas density. So it 
is naturally interesting to studying the potentials with variate density profile. 
A potentials with periodically variate profile were studied independently in 
[9], [10] in connection with a possibility of the parametric resonance in two 
neutrino system. In [9] an approximate solution for sinusoidal matter den- 
sity profile was found. In [10] the exact analytic solution for the periodic 
step-function density profile was obtained. In [TTJ an analytic approximate 
solution for arbitrary electron density for small matter effects on the neutrino 
oscillations was wound. 

In the present paper we have constructed toy model of electron density 
fluctuations. An approximate solution for two flavor neutrino oscillations 
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probabilities was found in periodic square potential profile with suggestion 
that the length of each square is smaller than the oscillation length. 

Neutrino oscillation in the periodic square potential 

In this section, we consider the mixing of two (active) neutrinos. Vf = 
Uv-mass) where Uf = (is e , u a ) are flavor states, v a is a linear combination of 
and v T , Vmass = iy\i lJ 2) T are mass states with masses m,\ and m 2 respectively. 

tt — ( cos ^ s ^ n ^ ^ (-\\ 
= { -sm6 cos9 ) [ } 

is a lepton mixing matrix in weak charged interactions. The equation for 
neutrino oscillations in matter is given by Schrodinger equation 
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Here V = y/2GpN e , N e is electron number density, Gf is Fermi constant, E 
is the energy of a neutrino, Am 2 = m 2 , — mf, x- the length of the neutrino 
beam path. 

We will assume that the V is of the form periodic rectangularly shaped 
potentials, as shown in Fig.l. 
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Figure 1: The dependence of the electronic potential V on distance L. Vb - 
is the averaged potential, AVis the amplitude of the additional rectangularly 
shaped potential and Ax is the half-period of the potential. 
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The formal solution of eq. (j2J), that is the oscillation matrix for the flavor 
states from initial point Xq to the final point Xf, can be written as 

S x ^ Xf = Te^ H ^ d \ (3) 

where T means the chronological ordering. 

Thus in order to get the oscillation matrix it is necessary to multiply 
consistently solutions on each period. Thus we get the oscillation matrix in 
the form of 

S= (e iC ) n = e inC , (4) 
where C is a hermitian matrix 

C = X -ln (y H+A V H ~ Aa; ) . (5) 

Here 

rr± tt( 0\ uT ^(Vo±AVO\ x f - x 

= V o o) u + [ o o) andAx = —^T (6) 

is the half period of the potential , n is the number of periods. 
The S matrix can be rewritten in the following form 

S = e™ c = e in ° tcdia3 ° = 0^e mCdia9 = OU ^ " . L ag | O, (7) 

i e m ° 22 I 

where C dia9 is diagonalized C matrix, C^x 9 and C^ 9 are its diagonal ele- 
ments, O is the matrix that diagonalize matrix C. 

0=( e 7 l)( C ° SIP V (8) 

y (J e 7 J \ — sin ip cos tp I 

The phase 7 is unphysical and does not enter into the transition probabilities. 
Solving these equation we find that 

• 2 4|Ci 2 | 2 , . 

sin zip = 2 2' l") 

(C11 — C22) + 4 I C12 1 



C diag_ c diag = ^ Cll - C 22 f + A\C 12 \ 2 . (10) 
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Now we can write the probability of v e — u a transition: 



p 



sin 2<^ sin n 



U 22 — U ll 



(11) 



We assume that the period of the potential, 2Ax, is smaller than the 
neutrino oscillation length in mater with average potential Vq and the am- 
plitude of the additional rectangularly shaped potential AV is smaller than 
(Ax) -1 . In that case we may perform series expansion for (jSJ) (according to 
the Hausdorf formula) 



C = C« + C< 2 > + + + o(Ax 5 ), 
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(13) 
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The hermitian matrix C has following elements 
Am 2 



C22 — C11 
Re C12 = 

Im C12 = 

Here 



E 



■cos29 m Ax + 0(Ax b 



~ sin20 m Ax- — ^sin26> m AV 2 Ax 3 + 0(Ax 5 ), 
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Am| 
4E 



12E 

Am 2 

sin 29 m AVAx 2 + ^ sin 29 m AV 

48E 

sin 29 



Am 2 
~2E~ 



AV 2 



(14) 

Ax 4 + 0(Ax 5 ). 



sin 29 r , 



(cos£-f^) +sin 2 # 



Ami 



Am \ cos 9 



2EV n " 2 



Am 2 



+ sin 2 9 



(15) 



(16) 



are the mixing angle and the effective mass square difference in the medium 
with the potential Vq. 
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Inserting (1141) into and (11 II) we obtain an approximate solution for 
the v e — > v a transition probability: 

(Am 2 \ 
—gi( Xf - Xo )(l + (] 2 )\, (17) 

where 

A = --^AV 2 Ax 2 cos 2 29 m + O (Ax 4 ) , (18) 
(3 2 = -^AV 2 Ax 2 sin 2 29 m + O (Ax 4 ) . (19) 

As it follows from eqs. (jTTj) - (jl~9~I) in the lowest approximation the v e — > v a 
transition probability is given by the MSW equation - the oscillations in the 
matter with the average potential Vo- 

In the first approximation, by the length of the period of the potential, 
the oscillation amplitude gets an additional factor ( JTBl and the oscillation 
phase shifts according to eqs. f|T7|) and ffl"9l) . Thus in that approximation only 
AV 2 'Ax 2 enters in the corrections. According to (fT4"j) . in the next approxi- 
mation the corrections (to the amplitude and to the phase of the oscillation 
probability) have forms (AV^^Ax 2 ) 2 and (AVAx) 4 . 

Conclusion 



An approximate analytic solution for two flavor neutrino oscillations prob- 
abilities was found in the periodic square potential profile with suggestion 
that the length of each square is smaller than the oscillation length in the 
matter with averaged potential and is smaller than the inverse amplitude 
((AV^)^ 1 ) of the additional rectangularly shaped potential. 
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